APCTP 97-05 
PM 97-05 
March 1997 



Higgs Boson Decays into Light Gravitinos 



Abdelhak Djouadi 1 and Manuel Drees 2 

1 Laboratoire de Physique Mathematique et Theorique, UPRES-A 5032, 
Universite de Montpellier II, F-34095 Montpellier Cedex 5, France. 

2 APCTP, 207-43 Cheongryangri-dong, Tongdaemun-gu, 
Seoul 130-012, Korea 



Abstract 

We study decays of the neutral and charged Higgs bosons of the minimal super- 
symmetric standard model into a very light gravitino and a neutralino or chargino. 
Present experimental constraints already imply that the lighter scalar Higgs boson 
can only have a very small branching ratio into this mode. However, if the gravitino 
mass is below ~ 0.5 eV, the heavier neutral and charged Higgs boson decays into a 
gravitino can be important or even dominant. We discuss the resulting signature for 
associate production of the heavy CP-even and the CP-odd states at future e + e~ 
colliders, and comment on Higgs boson production at hadron colliders. 



1. Introduction 



Recently models with a very light gravitino G [[|], mg < 10 -3 eV, have attracted some 
attention 0, [| ||. This interest was originally triggered by the resurgence of models 
II where supersymmetry breaking is transmitted by gauge interactions to the "visible 
sector" containing the usual gauge, Higgs and matter superfields, although these models 
tend to predict mg > 1 eV. Certain supergravity models can also naturally accommodate 
a very light gravitino 

Another motivation comes from CDF's observation of a single event with an e + e~ 
pair, two hard photons, and missing transverse momentum. The probability for such 
an event to come from Standard Model (SM) processes is about 10~ 3 , but it can quite 
easily be explained || in models with a light gravitino. It should be admitted that the 
absence of additional events with two hard photons and missing transverse momentum 
H is beginning to cast some doubt on this explanation. Nevertheless models with a 
very light gravitino remain an interesting field of study in their own right [p. 

In this Letter we explore the impact that a very light gravitino could have on the phe- 
nomenology of the Higgs sector of the minimal supersymmetric standard model (MSSM). 
Higgs bosons do not couple to pairs of gravitinos, but they do couple to a gravitino and a 
neutralino or a chargino. The corresponding branching ratios are computed in Section 2. 
Since in this scenario even the lightest neutralino is unstable, H — > G\ decays give rise to 
final states containing a hard photon, missing energy, and additional leptons and/or jets 
if X Xi- Higgs search strategies may have to be modified to include such final states. 
One example, the associate production of the heavier CP-even and CP-odd neutral Higgs 
bosons at a future e + e~ collider, is discussed in Section 3. Finally, Section 4 contains a 
brief summary and some conclusions. 



2. Decay branching ratios 

As first pointed out by Fayet [|l|, the couplings of the "longitudinal" (spin 1/2) components 
of the gravitino to ordinary matter are enhanced by the inverse of the gravitino mass; if 
is sufficiently small, this can compensate the suppression by the inverse Planck mass 
2.4 ■ 10 18 GeV that appears in all gravitational interactions. In fact, a longitudinal 
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nothing but the "Goldstino" that signals the spontaneous breaking 
of global supersymmetry (SUSY), and whose coupling are inversely proportional to the 
SUSY breaking scale A 2 ~ mgMp. Since Goldstino couplings contain momenta of the 
external particles, partial widths for decays into final states containing (longitudinal) 
gravitinos depend very strongly on the mass of the decaying particle. 



Specifically, for the case at hand, one has for the partial decay widths 
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where <3? stands for a neutral or charged Higgs boson and x% is one °f t ne f° ur neutralinos 
or two charginos of the MSSM. The coupling factors are given by 0: 



K ih = |Aj 3 sina; — Aj 4 cosa| 
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«iA = |iVi3 sin/5 + A" i4 cos/3 1' 
KiH± = |V^ 2 | 2 cos 2 /? + I fXj 2 1 2 sin 2 /?. (2) 



Following the notation of Ref . [12] , we have introduced the neutralino mixing matrix N, 
the two chargino mixing matrices U and V, the mixing angle a in the neutral CP-even 
Higgs sector, and the ratio of Higgs vacuum expectation values tan/3. The structure of 
eqs. (H) is easily understood, since gravitinos only couple to members of the same super- 
multiplet (in current basis). Each term is therefore the product of a higgsino component 
of a chargino or neutralino and the component of the corresponding Higgs current eigen- 
state in the relevant Higgs mass eigenstate. 

It would appear from eq. (H) than the partial widths for Higgs to gravitino decays could 
be made arbitrarily large by making mg very small, if m$ > m^. However, as stated 
earlier a very small mg corresponds to a small SUSY breaking scale A. As argued in 
Ref. 0, present lower bounds on sparticle masses imply that A should be several hundred 
GeV at least, which corresponds to mg > (a few) 10~ 5 eV. We therefore chose mg = 10~ 4 
eV in most of numerical examples. This corresponds to A = 650 GeV, which is probably 
quite close to its lower bound in realistic models. Note that A 2 is the vev of the largest 
F- or D-term in the model, which is usually (much) larger than the squared masses of 
sparticles in the visible sector. 

Furthermore, if mg < 1 eV, even the lightest sparticle in the visible sector will decay 
inside the detector |TJ, |2j . Here we make the usual assumption that the lightest neutralino 
Xi is this next-to-lightest sparticle, the LSP being of course the gravitino. Associate X\G 
production and XiXi P a i r production at e + e~ colliders then lead to final states with one 
or two isolated photons and missing energy-momentum. No significant excess over SM 
expectations has yet been observed in these channels. Following Refs. || |] we interpret 
this as implying 

(x(e + e- -> Gxl) < 0.1 pb (s/s = 91 GeV, | cos# 7 | < 0.7); (3a) 
a(e + e- -> XiXi) < 0.2 pb (y/s = 162 GeV). (3b) 



Expressions for these cross sections can be found in Refs. || and [|13|], respectively. In 
particular the LEP1 constraint ( |3a|) is very restrictive ||, if m^o < Mz- 

Further constraints on the relevant parameter space can be derived from the negative 
outcome of chargino searches. Note that in our scenario chargino pair production always 
leads to large amounts of visible energy, since Xi decays inside the detector. A small 
xf — Xi mass difference is thus no longer problematic for these searches. We have therefore 
interpreted recent chargino searches at LEP as implying 

m-± > 80 GeV. (4) 

In fact, in this scenario existing Tevatron data most likely imply a significantly stronger 
bound on the chargino mass. The authors of Ref. estimate that the present data 
samples of the CDF and DO collaborations should contain several additional events with 
two hard photons and substantial missing transverse momentum if m-± < 120 to 130 
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GeV. No such events have been observed ||. Although no "official" experimental bound 
on m~± has yet been published for the light gravitino scenario, it seems very likely that 
combinations of parameters giving m-± < 120 GeV are already excluded. We will there- 



fore indicate the consequences of replacing the lower bound 
constraint. 



by this more restrictive 
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Figure 1: Contour of constant BR(h — > Gx°) = 5% in the (/j,, M%) plane, for rrif L = 
rrif R = 1 TeV, A t = v6 TeV, tan/5 = 50, ttia = 500 GeV and mg = 10~ 4 eV. The region 
to the right and below the solid line is excluded by the experimental bounds (0) and (f|), 
where we have assumed = mg R = 300 GeV. In the region to the right and below 
the dotted contour the mass of the lighter chargino is below 120 GeV; this region is most 
likely excluded by present Tevatron data, as discussed in the text. 



We are now ready to present some numerical results. We begin with a discussion of 
the decays of the lighter neutral Higgs scalar h. Eq. (JJ) shows that T(h — > Gx°) will 
be maximal if is maximal. In Fig. 1 we have therefore chosen tan/3 = 50, a heavy 
pseudoscalar Higgs boson = 500 GeV, large degenerate third generation squark masses 
( m i L = m t R = 1 TeV), and sizable ti, — £r mixing (A t = \/6 TeV). This maximizes [OJ 
one-loop corrections to m^, which have been computed using the formalism of Ref. |I5| ; we 
obtain rrih = 139 GeV for our choice mt (pole) = 176 GeV.Q Here and in the subsequent 

*We include the leading 2-loop QCD corrections by using the running top mass at scale 
Q = ^Jm t mi in the expressions for the corrected Higgs boson masses. However, we neglect electroweak 
corrections JlTj , which could reduce mh by several GeV. 
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figures we assume that the U(l)y and SU{2) gaugino masses satisfy the "unification 
condition" M\ = 5/3 tan 2 8wM 2 — 0.5M 2 . Relaxing this assumption will not change our 
results significantly, since we are interested in the higgsino-like neutralino and chargino 
states, see eqs. (0). Further, we assume mg £ = rrie R = 300 GeV. In this case the constraint 
( j3b"D is always weaker than the combination of the constraints fl3a[) and (H). Our calculation 



of the total decay widths of the MSSM Higgs bosons includes decays |L8| into massive // 
pairs (/ = c, r, b, t), loop decays into two gluons or two photons and, where appropriate, 
decays into VV final states (V = W ± or Z, one of which may be off-shell) as well as 
decays of the heavy Higgs bosons into lighter ones (H — > hh, A — > Zh). 

The region to the right and below the solid line in Fig. 1 is excluded by present 
constraints; in the region of small the constraint (Q) is the relevant one, while for 
larger (gaugino-like xT) the condition (|3a|) is most restrictive. This still allows a 
narrow strip of parameter space where BR{h Gx°) exceeds 5%, as indicated by the 
long-dashed contour. However, this region can be excluded in the likely case that existing 
Tevatron data already impose the limit m^± > 120 GeV; this bound is violated to the 
right and below the dotted contour. We therefore conclude that, even if a light gravitino 
exists, it cannot change the h decay patterns significantly. 
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Figure 2: Contours of constant BR(A — > Gx°) m the (//, M 2 ) plane, for = 200 GeV 
and tan/3 = 1.5. The sfermion and gravitino mass parameters are as in Fig. 1. The region 
to the right and below the solid line is excluded by the constraints (f|) and while to 
the right and below the dotted contour the lighter chargino mass is below 120 GeV. 
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However, the situation could be quite different for the heavier Higgs bosons. As an 
example, in Fig. 2 we show contours of constant BR(A — ► Gx°) for tan/5 = 1.5 and 
rriA = 200 GeV; the values of the other parameters are as in Fig. 1. We have chosen a 
small value of tan/3 here in order to minimize the partial widths for A — > bb and A — > t + t~ 
decays. Since m 4 ~ and A t are large, ~ 100 GeV so that A — > Z/i decays are phase 
space suppressed. We see that in this case the gravitino mode can account for 40% of all 
A decays even if we require m-± > 120 GeV (dotted contour). This is true even though 
we chose a rather moderate value of in this figure, so that AH pairs can be produced 
at the proposed next linear e + e~ collider operating at y/s = 500 GeV. 

Before we discuss AH production in more detail, we explore how the results of Fig. 2 
depend on the values of various parameters. Increasing tan/3 tends to reduce the impor- 
tance of A decays into sparticles, at least for < 2m t , since the partial widths into bb 
and t + t - pairs scale oc tan 2 (3. Choosing the opposite sign of fi also reduces the maximal 
branching ratio for A — > G\° decays. The partial widths for these decays do not change 
much when the sign of fi is flipped, but the constraint (^) is more restrictive for fi > 0; 
in addition, the partial widths into pairs of neutralinos or charginos are enhanced. As a 
result, we find BR(A — > Gx°) < 25% for /i > 0, keeping all other parameters as in Fig. 2; 
this is reduced to 15% if we require m^± > 120 GeV. 




Figure 3: Branching ratios of the heavy Higgs bosons of the MSSM as a function of 
the mass of the pseudoscalar Higgs boson, for M 2 = 300 GeV, \i = —150 GeV, tan/3 = 
2, mg = 10 -4 eV and the same sfermion mass parameters as in Fig. 1. Here / stands for 
any SM quark or lepton, x° f° r an y of the four neutralinos of the MSSM, and x f° r an y 
of the two charginos or four neutralinos. 
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As discussed earlier, we probably cannot enhance the partial widths into gravitinos 
by reducing mg even further, since this would give too small a SUSY breaking scale A. 
However, we saw in eq. (|J) that these partial widths also rise very quickly with increasing 
Higgs boson mass. This is illustrated in Fig. 3, which shows some branching ratios of the 
heavy Higgs bosons of the MSSM as a function of tua- Here we have chosen M2 = 300 
GeV, n = —150 GeV and tan/3 = 2, corresponding to m^o = 136 GeV and m-± = 158 
GeV. The parameters of the squark mass matrices are as in Figs. 1 and 2. We see that, 
in spite of the m\ factor, the branching ratio for A — > Gy° decays at first grows rather 

slowly from the threshold at rn A = ny, due to the kinematical factor (l - ,„ 
eq. fll]). This is to be contrasted with the rapid increase of the branching ratio into heavy 
SM fermions just beyond the ti threshold (short dashed curve); here the threshold factor 

1 /9 

is (1 — Amf/m\) . Finally, the long dashed curve shows the branching ratio for A decays 
into pairs of neutralinos or charginos, which exhibits several thresholds. Together with the 
onset of A — > ti decays, these supersymmetric decay modes at first reduce BR(A — > Gx°) 
as rriA is increased beyond 320 GeV. These XX modes are also responsible for the reduced 
importance of the gravitino modes at small \fi\ shown in Fig. 2. However, eventually the 
m\ factor wins out, causing the solid curve to rise again as ttla is increased beyond 400 
GeV. 

So far we have only discussed decays of h and A. For large tua, tt? a ^ f- 2 , the partial 
width for H —>■ Gx° decays is very similar to that for A — > Gx° decays. However, if tan/? 
is not large, H can also have a sizable branching ratio into hh pairs JT8|. As shown by 



the dot-dashed curve in Fig. 3, this reduces the branching ratio for decays of the heavy 
neutral scalar into gravitinos. In fact, for the parameters of Fig. 2 this mode completely 
dominates H decays, leading to BR(H — > Gx°) < 2% in this case. On the other hand, 
for large the charged Higgs boson has an even slightly larger branching ratio into 
the gravitino mode than the pseudoscalar (dotted curve). However, at small tua, the 
branching ratio for H + — > G\ + decays is small, partly because the lightest chargino is 
heavier than the lightest neutralino, and partly because the threshold for H + — > tb decays 
is at much lower values of rriA than that for A — > ti decays. It should be emphasized that 
decays into light gravitinos could dominate the decays of all three heavy Higgs bosons of 
the MSSM, if uia > 700 GeV. Note also that for our choice mg = 10~ 4 eV, the total 
decay width of these Higgs bosons is around 100 GeV for ttia = 1 TeV. 

Finally, while it is unlikely that mg is significantly smaller than 10~ 4 eV, it could 
certainly be larger. We find that for mg = 5 ■ 10~ 4 eV, BR(A —>■ Gx°) < 15% for rriA < 1 
TeV. Clearly this is about the largest value of mg where this branching ratio can be 
sizable, if rriA lies in the range indicated by naturalness arguments. 



3. Collider Signals 

The signature for A — > Gx° decays depends on which neutralino states are produced 
in association with the light gravitino. Recall that these decays proceed through the 
higgsino components of the neutralinos. In most cases two of the four neutralinos of the 
MSSM are higgsino-like and the other two are gaugino-like, although occasionally two 
states are strongly mixed. If the higgsino-dominated states are the heavier neutralinos 
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X3 4> A — > Gx° decays will lead to the subsequent cascade decays X34 — ► where 
/ stands for any kinematically accessible quark or lepton; these decays proceed through 
(virtual or real) h, Z and / exchange. X34 — > Xiff decays may also occur. This 
can lead to quite complicated final states. However, in this scenario some decays of the 
type A — > Xi,2X3,4 are probably also allowed. Since they combine a higgsino-like and a 
gaugino-like state, the relevant couplings are large [fL2]l . In this case A — > Gx° decays 
are therefore only important if itia is quite large, probably beyond the reach of next 
generation colliders. 

We therefore focus on the scenario where the higgsino-like states are lighter than 
the gaugino-like states. We saw in Fig. 2 that in such a scenario the branching ratio 
for A — > Gx° decays can be substantial already for = 200 GeV. Note that the 
mass splitting between the higgsino-like states is usually quite small. A — > Gx® decays 
therefore differ from A — > G\\ decays only by a rather soft // pair. With the possible 
exception f = fi, such soft pairs are probably not detectable at all at hadron colliders like 
the planned LHC. At these colliders the most important Higgs production mechanism is 
gluon fusion, gg — > h, H, A. The signature for A — > G\° decays would then be a single 
hard photon from — > 7G and missing transverse momentum. Note that there is no 
Jacobian peak in the Pt(j) distribution, since the photon is only produced in a secondary 
decay. Unfortunately, this final state will probably suffer from considerable backgrounds^. 
Recall also that in this scenario all SUSY processes give rise to two hard photons and 



missing @, |10[- They can therefore also contribute to the background if one of the 
photons is lost.[] A detailed study will be necessary to decide whether vetoing the presence 
of additional jets and/or leptons will be sufficient to overcome this potentially very large 
background. 

As usual, things should be much "cleaner" at e + e~ colliders, especially if A is light 
enough to be produced at a 500 GeV collider. The main A production mechanism at 
e + e~ colliders is associate AH production. We saw in Sec. 2 that if rriA < 220 GeV, so 
that the cross section for AH production is sizable at a 500 GeV collider, and if tan/3 
is relatively small, so that BR(A — > Gx°) can be substantial, then BR(H — > Gx°) is 
very small, largely due to H — > hh decays. AH production and subsequent A — > Gx° 
decay would then most often lead to a final state with two bb pairs from the decays of 
the two light Higgs scalars, a hard isolated photon, and perhaps a soft // pair from 
X2 — > Xiff decays. There should be little background in this channel, especially if one 
requires several b— tags. Furthermore, events of this type should allow to determine the 
masses of all involved (s)particles by studying kinematic distributions. Clearly rrih can be 
obtained from bb invariant mass distributions. Combinatorial backgrounds can be reduced 
by choosing the pairing that gives the most similar values for the two bb invariant masses. 
Of course, rrih can also be obtained completely independently from e + e~ — ► Zh events in 



TFor instance as irreducible backgrounds, one has "/Z production and subsequent Z — > vv decay; W*y 
final states followed by W — > li>\ (I = e, fi, t) decays can also contribute to the background if the lepton 
is not detected. 

*For m~ — eV, q — > Gq decays will dominate even over q — > gq decays if rriq > 800 GeV 

similarly, gluinos will dominantly decay into a gluon and a gravitino if m~ < 10~ 3 eV • (rriq/1 TeV) 2 . 
These decays lead to final states with rather low jet multiplicity, increasing the probability that the jets 
escape detection. 
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this scenario; the value of derived from the kind of AH events we are studying here 
would therefore be a cross-check on the interpretation of the observed events. 

In principle one can get m# from the invariant mass of the 4b system, and itla as the 
mass recoiling against this system. However, the b — > c — > s cascade decays of the four 
b (anti) quarks are likely to produce some neutrinos, which will smear out these invariant 
mass distributions. Once m/, is known, one might be able to at least partly correct for 
this by using constraint fits. Alternatively, one can obtain tua and mn, as well as m 
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from the photon energy spectrum. Two examples of this spectrum are shown in Fig. 4, 
for m A = 200 GeV, m H = 220 GeV, y/s = 500 GeV and m^> = 100 GeV (solid) and 150 

GeV (dashed). Here we have treated G\\ final states the same as Gxi final states; the 
additional soft // pair does not distort the E 1 spectrum significantly. 
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Figure 4: The photon energy spectrum in the lab frame for events of the type e + e _ — > HA 
with subsequent A — > Gxi decay and Xi ~^ G^y. We have taken itlh = 220 GeV, = 200 
GeV, y/s = 500 GeV, and two different values of the Xi mass, as indicated. The acceptance 
cut | cos# 7 | < 0.9 has been applied. 



We see that the photons are always quite hard; an acceptance cut like _E 7 > 10 GeV 
would therefore not reduce the signal at all. In Fig. 4 we have required |cos# 7 | < 0.9, 
where # 7 is the angle between the outgoing photon and the e~ beam direction in the 
lab frame. This reduces the accepted cross section by some 8%. Recall that da(e + e~ — > 
AH)/dcos0A oc sin 2 ^, where 9a is the angle between the pseudoscalar Higgs boson and 
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the e beam; the distribution in cos# 7 is therefore mildly peaked at cos# 7 = 0. 

Of course, this angular acceptance cut has no bearing on the endpoints of the photon 
energy spectrum, which are at 



where -Emax(Xi) ^ s the maximal energy of the neutralino in the lab frame and Pmax(Xi) 

-^max(Xi) — m \ is the absolute value of the corresponding 3-momentum. Since A — > Gx° 
decays are isotropic, the distribution in E(xi) is A^, with endpoints 

£ m ax,min(X?) = — \E { A)(X°i) ± kE [A) (G)] . (6) 

Here, 

m \ + ml ~ m ^ — mi 

are the neutralino and primary gravitino energy in the A rest frame; Ea = (s + m^ — 



m 2 H ) / (2yfs) is the energy of A in the lab frame and k = J E\ — m? A is the absolute value 
of the 3-momentum of A (or H) in the lab frame. 

The determination of the endpoints of the photon energy spectrum already gives two 
constraints on the three unknowns tua, itlh and m^o. In addition, Fig. 4 shows that the 
photon energy spectrum is almost flat over some range of energies; indeed, it would be 
completely flat had we not imposed the acceptance cut | cos# 7 | < 0.9. This is due to 
the fact that both the A — > Gx° decay and the x\ — > G7 decay are isotropic in their 
respective rest frames. In the absence of acceptance cuts this leads to a flat spectrum 
between the points 

E^^[E min {x 1 )±p min {x° 1 )}, (8) 

where E min (xi) has been given in eq. ([]). The solid curve in Fig. 4 shows that this plateau 
can be very narrow. This happens for combinations of parameters leading to a small value 
of Pmin(x?); hi the given case this quantity amounts to just 6 GeV. In more generic cases, 
one should be able to determine both E^ and E~ with good precision, given sufficient 
statistics. This would give two additional, independent constraints on three unknown 
masses. 

Another possibly useful kinematic distribution is that in the invariant mass of the two 
gravitinos produced in the A — > Gx° — > GG7 decay chain. Since the decays are isotropic, 



one simply finds da/dM^Q oc Mgg with a sharp cut-off at M~~? = — m|. However, 

the GG invariant mass will be equal to the missing mass only if there are no additional 
sources of missing energy in the event (neutrinos from the decay of b quarks, or initial 
state radiation that goes down one of the beam pipes). 
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4. Summary and Conclusions 

In this note we have studied the decays of MSSM Higgs bosons into light gravitinos. We 
found that present experimental bounds already exclude the possibility that such decays 
are relevant for the light neutral Higgs scalar, but they could be quite important, or even 
dominant, for the heavier Higgs bosons if mg < 0.5 meV. In the simplest case this gives 
rise to final states with a hard isolated photon and missing (transverse) energy. It is 
doubtful whether this signal is detectable at the LHC, but associate AH production at 
e + e _ colliders will have very little background if only A decays into a gravitino and a 
neutralino while H decays into a fully reconstructable final state. We argued that such 
a scenario is in fact likely if A is sufficiently light to be produced at a 500 GeV collider. 
We also showed that in this favorable situation one can in principle reconstruct m#, 
and m^o from the energy spectrum of the photon. 

In practice it will probably be better to reconstruct the neutralino mass directly from 
neutralino pair events. Since the neutralino must have a substantial higgsino component 
if the branching ratio for A —>■ Gx° decays is sizable, the cross section for neutralino pair 
production at the same e + e~ collider will be quite substantial - in fact, significantly larger 
than the AH production cross section. This also implies that one need not worry about 
Higgs decays into gravitinos if one does not find sizable signals for neutralino and chargino 
pair production at the same energy. In this respect Higgs decays into light gravitinos are 
quite similar to Higgs decays into pairs of neutralinos or charginos. At present, however, 
all these possibilities must be kept in mind when deriving search strategies for Higgs 
bosons at future colliders. 

Acknowledgements: M.D. thanks the members of the LPM at the Universite de Mont- 
pellier for their hospitality during a visit where this project was initiated. 
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